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Abstract — This work presents joint interference suppression 
and power allocation algorithms for DS-CDMA networks with 
multiple hops and decode-and-forward (DF) protocols. A scheme 
for joint allocation of power levels across the relays subject to 
group-based power constraints and the design of linear receivers 
for interference suppression is proposed. A constrained minimum 
mean-squared error (MMSE) design for the receive filters and the 
power allocation vectors is devised along with an MMSE channel 
estimator. In order to solve the proposed optimization efficiently, 
a method to form an effective group of users and an alternating 
optimization strategy are devised with recursive alternating least 
squares (RALS) algorithms for estimating the parameters of the 
receiver, the power allocation and the channels. Simulations show 
that the proposed algorithms obtain significant gains in capacity 
and performance over existing schemes. 

Index Terms — DS-CDMA, cooperative systems, optimization 
methods, adaptive algorithms, resource allocation. 



I. Introduction 

Multiple-antenna wireless communication systems can ex- 
ploit the spatial diversity in wireless channels, mitigating 
the effects of fading and enhancing their performance and 
capacity. Due to the size and cost of mobile terminals, it is 
considered impractical to equip them with multiple antennas. 
However, spatial diversity gains can be obtained when single- 
antenna terminals establish a distributed antenna array via 
cooperation fl1-||3l. This allows a significant reduction on 
the transmitted power for an equivalent performance. In a 
cooperative system, terminals or users relay signals to each 
other in order to propagate redundant copies of the same 
signals to the destination user or terminal. To this end, the 
designer must use a cooperation protocol such as amplify- 
and-forward (AF) [3|, decode-and-forward (DF) fl3], |@] and 
compress-and-forward (CF) [5|. 

Recent contributions in the area of cooperative and multihop 
communications have considered the problem of resource 
allocation 0, (7J. Prior work on cooperative multiuser DS- 
CDMA networks has focused on the assessment of the impact 
of multiple access interference (MAI) and intersymbol inter- 
ference (ISI), the problem of partner selection [4], [8|, the bit 
error ratio (BER) and outage performance analysis [9], and 
training-based joint power allocation and interference mitiga- 
tion strategies [10|, [12|. However, these strategies require a 
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higher computational cost to implement the power allocation 
and a significant amount of signalling, decreasing the spectral 
efficiency of cooperative networks. This problem is central to 
ad-hoc and sensor networks [13] that utilize spread spectrum 
systems and require multiple hops to communicate with nodes 
that are far from the source node. 

In this work, joint interference suppression and power 
allocation algorithms for DS-CDMA networks with multiple 
hops and DF protocols are proposed. A scheme that jointly 
considers the power allocation across the relays subject to 
group-based power constraints and the design of linear re- 
ceivers for interference suppression is proposed. The idea of 
a group-based power allocation constraint is shown to yield 
close to optimal performance, while keeping the signalling 
and complexity requirements low. A constrained minimum 
mean-squared error (MMSE) design for the receive filters 
and the power allocation vectors is developed along with an 
MMSE channel estimator for the cooperative system under 
consideration. The linear MMSE receiver design is adopted 
due to its mathematical tractability and good performance. 
However, the incorporation of more sophisticated detection 
strategies including interference cancellation with iterative de- 
coding lfT4l - lfT6l and advanced parameter estimation methods 
ifTTll - ETI are also possible. In order to solve the proposed 
optimization problem efficiently, a method to form an effective 
group of users and an alternating optimization strategy are 
presented with recursive alternating least squares (RALS) 
algorithms for estimating the parameters of the receiver, the 
power allocation and the channels. 

The paper is organized as follows. Section II describes a 
cooperative DS-CDMA system model with multiple hops. Sec- 
tion III formulates the problem, details the constrained MMSE 
design of the receive filters and the power allocation vectors 
subject to a group-based power allocation constraint, and 
describes an MMSE channel estimator. Section IV presents an 
algorithm to form the group and the alternating optimization 
strategy along with RLS-type algorithms for estimating the 
parameters of the receiver, the power allocation and the 
channels. Section V presents and discusses the simulation 
results and Section VI draws the conclusions of this work. 

II, Cooperative DS-CDMA Network Model 

Consider a synchronous DS-CDMA network with multipath 
channels, QPSK modulation, K users, N chips per symbol 



and L as the maximum number of propagation paths for 
each link. The network is equipped with a DF protocol that 
allows communication in multiple hops using n r fixed relays 
in a repetitive fashion. We assume that the source node or 
terminal transmits data organized in packets with P symbols, 
there is enough control data to coordinate transmissions and 
cooperation, and the linear receivers at the relay and destina- 
tion terminals are synchronized with their desired signals . The 
received signals are filtered by a matched filter, sampled at chip 
rate and organized into M x 1 vectors r s d[mj), r sri [rrij] and 
r ri d[rrij\, which describe the signal received from the source 
to the destination, the source to the relays, and the relays to 
the destination, respectively, 

K 

r s d[mi] = ^2a^ d [mi]C k h s d,k[mi]bk[mi] + ri sd [mi] 
fc=i 

+ n sd [m 1 } 1 

K 

r srj [mx] = ^2 a sr 3 [rni]C k h srjtk [mi]b k [mi] + ri srj [m{\ 
fc=i 

+ n sr] [m 1 } 1 

K 

r rj d[mj] = ^a*. d [m j ]C k h rjd ,klm j }bk[m j ] + r} r . d [m j } 
fc=i 

+ n rjd [m J ] 1 

j = 1, . . . , n p , rrij = (j - 1)P + l,...,jP, 
i = l,...,P 

(1) 

where M = N + L — 1, P is the number of packet symbols, 
rtp = n r + 1 is the number of transmission phases or hops, n r 
is the number of relays, and rrij is the index of original and 
relayed signals. The vectors n sd \m\\, n sr . [mi] and n r .d[rrij\ 
are zero mean complex Gaussian vectors with variance er 2 
generated at the receivers of the destination and the relays 
from different links, and the vectors r) sd [mi\, r) sr [mi] and 
r) rd [ m j] represent the intersymbol interference (ISI). The 
amplitudes of the source to destination, source to relay and 
relay to destination links for user fc are denoted by aj d [mi], 
a k sr . [mi] and a k , d [mj], respectively. The quantities fe/Jmi] and 
bk[rrij] represent the original and reconstructed symbols by 
the DF protocol at the relays, respectively. The M x L matrix 
C k contains versions of the signature sequences of each user 
shifted down by one position at each column as described by 



c fe (l) 

! '■■ cfc(i) 

c k (N) ; 

■■■ c k (N) 



(2) 



where c k = [cfe(l), c k {2), . . . , c k (N)\ stands for the signa- 
ture sequence of user k, the Lxl channel vectors from source 
to destination, source to relay, and relay to destination are 
h s d,k[mi], h sri ,k[mi], h rjd> k[mj], respectively. By collecting 
the data vectors in (JTJ (including the links from relays to 
the destination) into a (n r + \)M x 1 received vector at the 



destination we obtain 

r s d[mi] 

r ri d[m 2 ] 



r d[m r . 



^K 



J2k=i a s d[mi]Ckh s d,k[mi]b k [i] 
J2k=i a ^ ld [ m 2\C kh ri d, k [ m 2]b 
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Rewriting the above signals in a compact form yields 

K 

" J2^k[i\A k [^khk[i\+v[i\ 



ri 



■n\i\ 



fc=i 



K 

E 

k=l 
K 



B k [i]A k [i]C k h k [i] + rj[i] + n[i] 



(4) 



J2 p k[i]B k [i]a k [i] + r][i] 



fc=i 



where 

e k = 



l)M 



(»< 



DL 



matrix 
shifted 



the (n r 

diag{Cfc . . . C k } contains copies of C k 
down by M positions for each group of L columns 
and zeros elsewhere. The [n r + 1)L x 1 vector h k [i] 
contains the channel gains of the links between the 
source, the relays and the destination, and p k [i] is the 
effective signature for user k. The (n r + 1) x (n r + 1) 



diagonal matrix B k [i] 



diag(M*] bl} a [i] 



■b r k d m 

contains the symbols transmitted from the source to the 
destination (6fe[i]) and the n r symbols transmitted from 
the relays to the destination (K 1 [«]... 61™ [«]) on the main 
diagonal, and the (n r + 1)M x (n r + 1)M diagonal matrix 



B k \i] = diag(6 fc [i]®I M b r k ld [i}®I 



M ■ 



lr n d\ 



\®I 



M 



where (££) denotes the Kronecker product and Im is an 
identity matrix with dimension M. The (n r + 1) x 1 power 
allocation vector a k [i] = [a k sd [mi] a k rid [m 2 ] . ■ ■ a k nrd [m, lp ]] T 
has the amplitudes of the links, the (n r + 1) x (n r + 1) 
diagonal matrix A k [i] is given by A k [i] — dia,g{a k [i]}, 
and the (n r + 1)M x (n r + 1)M diagonal matrix A k [i] = 
[a k sd [mi] <8> Im a k rid [m 2 ] ® I M . . . a k nrd [m np ] ® I M ] T . 
The (n r + 1)M x (n r + 1) matrix P k has copies of the 
effective signature p k [i] shifted down by M positions for 
each column and zeros elsewhere. The (n r + 1)M x 1 vector 
rj[i] represents the ISI terms and the (n r + 1)M x 1 vector 
n[i] has the noise components. 

III. Proposed MMSE Receiver Design, Power 
Allocation and Channel Estimation 

In this section, a joint receiver design and power allocation 
strategy is proposed using constrained linear MMSE estima- 
tion and group-based power constraints along with a linear 
MMSE channel estimator. To this end, the (n r + 1)M x 1 
received vector in (@) can be expressed as 



k=£S 



Ps\i\B s \i]as >k [i]+y ] P k \i]B k [i]a k [i]+r)[i]+n[i], 

(5) 

where S = {Si,S 2 , ■ ■ ■ ,Sq} denotes the group of G 
users to consider in the design. The (n r + 1)M x 



G(n r + 1) matrix P s = [P Sl P$ 2 ... Ps G ] con- 
tains the G effective signatures of the group of users. 
The G(n r + 1) x G(n r + 1) diagonal matrix Bs[i] = 
diag^W b r s \ d \i]...b r s - d \i] ... b So \i] ^[i-^W) 
contains the symbols transmitted from the sources to 
the destination and from the relays to the destination 
of the G users in the group on the main diagonal, 
the G(n r + 1) X 1 power allocation vector as,k [i] = 



L«fiH<\jl*J 



• a fLdW 



Sg r 



•) a sd 



*&l 



% r„ r dl 



Of 



the amplitudes of the links used by the G users in the group. 

A. Linear MMSE Receiver Design and Power Allocation 
Scheme with Group-Based Constraints 

The linear MMSE interference suppression 
for user k is performed by the receive filter 

w k [i] = [wk,i[i], ■•■, u>jfe,(n r +i)Ar[*]] with (n r + 1)M 
coefficients on the received data vector r[i] and yields 



ZkW 



wf\i]r\i], 



(6) 



where Zk [i] is an estimate of the symbols, which are processed 
by a slicer Q(-) that performs detection and obtains the desired 
symbol as b k [i] = Q(z k [i\). 

Let us now detail the linear MMSE-based design of the 
receivers for user k represented by Wk[i] and for the compu- 
tation of the G(n r + 1) x 1 power allocation vector os,*^]. 
This problem can be cast as 



[-■.op* n °P t ' 

[ w k J a S,k 



E[(\b k [i]-t 
subject to as k [i]as,k[i] = Pg 



arg mm 

w k [i\,a s ,k[i} 



M\ 



(7) 



The MMSE expressions for the receive filter Wk[i] and the 
power allocation vector 05 a, [i] can be obtained by employing 
the method of Lagrange multipliers with (0, which leads to 



C k = E[\b k [i\ - w^[{\(P s [i]Bs[i\a s , k [i 
+ J2 Ph[i]B h [i]a h [i]+r)[i] + n[i] 



(8) 



k ^s 
+ \ k (a s , k [i] 



Pa), 



where A^ is a Lagrange multiplier. An expression for 05 fc[i] 
is obtained by fixing tOfc[i], taking the gradient terms of the 
Lagrangian and equating them to zero, which yields 



a s ,k\i] = (Rs,k\i] + AfcJ) 1 p s ,k[i] 



(9) 



where the G(n r + 1) x G(n r + 1) covariance matrix 
Rs, k \i] = E[B%[i]P§\t\w k [i]wZ[i]Ps[i]B s [i]] and the 
vector p Stk \i] = E[b k [i]B§[i]P%[{\w k [i]] is a G{n r + 1) x 1 
cross-correlation vector. The Lagrange multiplier A& plays the 
role of a regularization term and has to be determined numer- 
ically due to the difficulty of evaluating its expression. Now 
fixing 05 fc[»], taking the gradient terms of the Lagrangian and 
equating them to zero leads to 



w k [i] =R 1 [i]Pk[i], 



(10) 



where the covariance matrix of the received vector is given 
by R[i] = E[r\i]r H \i]} and p k [i] = E[b* k \i]r\i}] is a 
(n r + \)M x 1 cross-correlation vector. The quantities R[i] 
and Pk[i] depend on the power allocation vector as , k [i\. The 
expressions in (O and ( TTOb do not have a closed-form solution 



as they have a dependence on each other. Moreover, the 
expressions also require the estimation of the channel vector 
h k [i]. Thus, it is necessary to iterate (0 and (fTOt with initial 
values to obtain a solution and to estimate the channel. The 
network has to convey the information from the group of 
users which is necessary to compute the group-based power 
allocation including the filter i«fc[i]. The expressions in (0 
and ( TTOb require matrix inversions with cubic complexity ( 
0(((n r + 1)M) 3 ) and 0((K(n r + l)) 3 ). 

B. Cooperative MMSE Channel Estimation 

In order to estimate the channel in the cooperative system 
under study, let us consider the transmitted signal for user k, 
x k [i] — B k [i]A k [i]C k h k [i] = Q k [i]h k [i], and the covariance 
matrix given by 

R=[r\i]r H \i]] 

K 

= £ Qk[i\E[hk[i]hf[i\]Qk[i\ + E[ V [i\r] H [z\] + o 2 I 
fc=i 

K 



= ^Q fc [z]P ht gf[z] + P, + a 2 / 



fc=i 



(11) 



A linear estimator of h k [i] applied to r[i] can be represented 
as h k [i] = T k r[i\. The linear MMSE channel estimation 
problem for the cooperative system under consideration is 
formulated as 



T k , opt = a,Tgmin E[\\h k \i] - h k \i}\\ 2 ] 

Tk 

= d xgminE[\\h k [i)-T k i r\i}\\ 2 }- 



(12) 



Computing the gradient terms of the argument and equating 
them to zero yields the MMSE solution 



fc,opt 



R x Pk, 



(13) 



where P k = E[r[i]h£[i]] ^ = Q k $E[h k [i]hj[[i]] = 
Q k [i]Ph k - Using the relation h k [i] = T k r[i], we obtain 

hk\i}=T% opt r[i\ = P k i R- 1 r{i\ 

K 

= PhM [*1(E Qf[i\ p huQk W + Pn + ^iy'rii], 
fc=i 

(14) 

The expressions in (fl4] > require matrix inversions with cubic 
complexity ( 0(((n r + l)M) 3 )), however, this matrix inversion 
is common to dTOb and needs to be computed only once for 
both expressions. In what follows, computationally efficient 
algorithms with quadratic complexity (0(((n r + 1)M) 2 )) 
based on an alternating optimization strategy will be detailed. 

IV. Proposed Adaptive Algorithms 

In this section, we develop adaptive RALS algorithms using 
a method to build the group of G users based on the power 
levels, and then we employ an alternating optimization strategy 
for efficiently estimating the parameters of the receive filters, 
the power allocation vectors and the channels. Despite the joint 
optimization that is associated with a non-convex problem, the 



proposed RALS algorithms have been extensively tested and 
have not presented problems with local minima. 

The first step in the proposed strategy is to build the group 
of G users that will be used for the power allocation and 
receive filter design. A RAKE receiver is employed to obtain 
zf AKE [i] = (C k h k [i]) H r[i\ = pf [i]r[i] and the group is 



formed according to 
compute the G largest 



RAKEr 



1,2,. 



(15) 
The design of the RAKE and the other tasks require channel 
estimation. The power allocation, receive filter design and 
channel estimation expressions given in (0, (TTOb and (fPfl i. 
respectively, are solved by replacing the expected values with 
time averages, and RLS-type algorithms with an alternating 
optimization strategy. In order to solve (fl4t efficiently, we 
develop a variant of the RLS algorithm that is described by 



H 



h k [i\ = p hh [i\Q%[i\R 



\i\r\i\ 



(16) 



where Q k 



Bk[i]Ak[i]Ck, the estimate of the inverse of 



the covariance matrix R [i] is computed with the matrix 
inversion lemma 



k[%\ = 



- l R[i- l]r[i] 



l + a -ir H [i]R[i-l]r[i]' 



R\. 



l R[i - 1] - a- l k[i)r H [i)R[i - 1] 



and 



-H, 



P hh [i] = aP hk [i - 1] + h k [i - l]h k [i - 1], 



(17) 
(18) 

(19) 



where a is a forgetting factor that should be close to but less 
than 1, The approach for allocating the power within a group 
is to drop the constraint, estimate the quantities of interest and 
then impose the constraint via a subsequent normalization. The 
group-based power allocation algorithm is computed by 

a s ,k[i\ =-Rs,fc[*]Ps,k[*] 

= fls, fc [i](aps,fc[* - 1] + h\i}v k \i}) (20) 
= as,k[i- 1] +£a[*]*:s,fc[*], 

where £ a [i] — bk[i] — ^s fc[* — l] t, k[*] is the a priori error, 
Vk[i] = B s [i]Pg [i] w k[i] is the input signal to the recursion 



k 



sM l \ 



l Rs,k[*- 1 H[*\ 



1 



-1„H 



v^[i]R s , k [i-l]v k [i\ 



(21) 



Rs,k[i]=a 1 Rs,k[i-l]-a 1 k s , k [i]v"[i]R s . k [i -1]. 

(22) 
The normalization 05^ [i] <— Pg o,s,k[i]/\\as,k\i\\\ is then 
performed to ensure the power constraint. The receive filter is 
computed by 



Wk[i] = w k [i - 1] + fc[i]^*[«], 



(23) 



where the a priori error is given by £[z] = b k [i] — w k [i — 
l]r[i] and k[i] is given by {TJ\ . The proposed scheme employs 
the algorithm in ( PT31 l to allocate the users in the group and 
the channel estimation approach of (I16t -(|19t, The alternating 
optimization strategy uses the recursions (l20l and (l23~t with 
1 or 2 iterations per symbol i. 



V. Simulations 

The bit error ratio (BER) performance of the proposed joint 
power allocation and interference suppression (JPAIS) scheme 
and RALS algorithms with group-based power constraints 
(GBC) is assessed. The JPAIS scheme and algorithms are 
compared with schemes without cooperation (NCIS) and with 
cooperation (CIS) [8| using an equal power allocation across 
the relays (the power allocation in the JPAIS scheme is 
disabled). A DS-CDMA network with randomly generated 
spreading codes and a processing gain N ~ 16 is considered. 
The fading channels are generated considering a random power 
delay profile with gains taken from a complex Gaussian 
variable with unit variance and mean zero, L = 5 paths 
spaced by one chip, and are normalized for unit power. The 
power constraint parameter PA, k is set for each user so that 
the designer can control the SNR (SNR = Pa^/o 2 ) and 
Pt = Pq + (K — G)PA, k , whereas it follows a log-normal 
distribution for the users with associated standard deviation 
equal to 3 dB. The DF cooperative protocol is adopted and 
all the relays and the destination terminal use either linear 
MMSE, which have full channel and noise variance knowl- 
edge, or adaptive receivers. The receivers are adjusted with 
the proposed RALS with 2 iterations for the JPAIS scheme, 
and with RLS algorithms for the NCIS and CIS schemes. 
We employ packets with 1500 QPSK symbols and average 
the curves over 1000 runs. For the adaptive receivers, we 
provide training sequences with 7V tr = 200 symbols placed 
at the preamble of the packets. After the training sequence, 
the adaptive receivers are switched to decision-directed mode. 

The first experiment depicted in Fig. Q] shows the BER 
performance of the proposed JPAIS scheme and algorithms 
against the NCIS and CIS schemes with n r = 2 relays. The 
JPAIS scheme is considered with the group-based power con- 
straints (JPAIS-GBC). All techniques employ MMSE or RLS- 
type algorithms for estimation of the channels, the receive 
filters and the power allocation for each user. The results show 
that as the group size G is increased the proposed JPAIS 
scheme and algorithms converge to approximately the same 
level of the cooperative JPAIS-MMSE scheme reported in [ 10|, 
which employs G = K for power allocation, and has full 
knowledge of the channel and the noise variance. 

The proposed JPAIS-GBC scheme is then compared with a 
non-cooperative approach (NCIS) and a cooperative scheme 
with equal power allocation (CIS) across the relays for 
n r =1,2 relays. The results shown in Fig. [2] illustrate the 
performance improvement achieved by the JPAIS scheme and 
algorithms, which significantly outperform the CIS and the 
NCIS techniques. As the number of relays is increased so 
is the performance, reflecting the exploitation of the spatial 
diversity. In the scenario studied, the proposed JPAIS-GBC 
with G = 3 can accommodate up to 3 more users as compared 
to the CIS scheme and double the capacity as compared with 
the NCIS for the same BER performance. The curves indicate 
that the GBC for power allocation with only a few users is 
able to attain a performance close to the JPAIS-GBC with 
G = K users, while requiring a lower complexity and less 
network signalling. A comprehensive study of the signalling 
requirements will be considered in a future work. 
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Fig. 1. BER performance versus number of symbols. Parameters: Ay = 
\ k = 0.025 (for MMSE schemes), a = 0.998, Rs\[i] = 0.01 J and 
-R -1 [i] =0.017. 
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Fig. 2. BER performance versus SNR and number of users for the optimal 
linear MMSE detectors. Parameters: a = 0.998, Kg fe [i] = 0.01/ and 
R _1 fil = 0.017. 



VI. Concluding Remarks 

This work has proposed the JPAIS scheme with group-based 
constraints (GBC) for cooperative DS-CDMA networks with 
multiple hops and the DF protocol. A constrained MMSE de- 
sign for the receive filters and the power allocation with GBC 
has been devised along with an MMSE channel estimator. We 
have proposed RALS algorithms for estimating the parameters 
of the channels, the receive filter and the power allocation. The 
results have shown that the JPAIS scheme with GBC and the 
RALS algorithms achieve significant gains in performance and 
capacity over existing schemes. 
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